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WATER  QUALITY  MONITORING  AT  THE  NORFOLK  DISPOSAL  SITE 

By 

Raymond  W.  Alden1,  Joseph  H.  Rule2,  Suzanne  S.  Jackman3,  and  Arthur  J.  Butt4 


INTRODUCTION 

The  continual  dredging  of  navigational  channels  in  major  seaports  is 
essential  to  maintain  shipping  access.  However,  the  disposal  of  potentially 
contaminated  dredged  materials  raises  environmental  issues.  Current  methods 
for  dredged  material  disposal  include:  landfill,  onshore,  and  open  ocean 
disposal.  Available  land  for  onshore  and  landfill  disposal  is  often  at  a 
premium  in  the  industrialized,  urban  seaport  and  may  present  a  variety  of 
social,  economic,  and  ecological  problems. 

Recently,  renewed  interest  has  been  generated  in  the  feasibility  of  the 
open  ocean  disposal  of  dredged  materials.  Currently,  dredged  materials  from 
the  Hampton  Roads  Harbor  area  are  being  disposed  at  the  Craney  Island  con> 
tainment  facility.  However,  Craney  Island  has  a  finite  capacity  in  its 
current  configuration.  An  open  ocean  disposal  site  designated  the  Norfolk 
Disposal  Site  (NDS)  is  being  considered  as  an  alternative  for  some  of  the 
Norfolk  Harbor  system  dredged  material.  The  present  ongoing  study  involves 
the  monitoring  of  baseline  water  quality  characteristics  at  the  NDS.  This 
study  was  undertaken  to  characterize  the  magnitude  of  natural  spatial- tempo¬ 
ral  variability  of  various  ecological  parameters  and  to  develop  a  series  of 

r'  f' 

multivariate  statistical  models  to  be  used  as  an  ‘'early  warning  system"  in 
historic  trend  assessment  studies.  This  report  presents  the  results  of  a 
three  year  baseline  water  quality  program  at  the  NDS.  ^ - 
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STUOY  AREA 


The  proposed  Norfolk  Disposal  Site  delineated  as  a  circle  with  a  4  nm 
(7.4  km)  radius  is  located  beyond  the  10  fathom  contour  line  approximately 
27  Km  east  of  Chesapeake  Bay  mouth  (Figure  1).  The  water  flow  pattern  in 
this  region  of  the  middle  Atlantic  Bight  is  variable  (Beardsley  and 
Boicourt,  1980;  Boicourt,  1981).  The  flow  pattern  on  the  continental  sr.elf 
is  typically  southward  (Bumpus  1973)  with  the  shallower  inner  shelf  domi¬ 
nated  by  wind  driven  forces  (Boicourt  and  Hacker,  1976;  Boicourt,  1981). 
Near  estuarine  influences,  the  flow  patterns  are  more  complex.  More  speci¬ 
fically,  the  flow  pattern  at  the  Chesapeake  Bay  mouth/conti nental  shelf 
interface  and  seaward  is  very  dynamic.  Circulation  at  the  interface  is  in 
response  to  the  synergistic  interactions  of  river  run-off,  vertical  decoup¬ 
ling  at  the  pycnocline,  wind  and  tidal  prism  patterns  (Boicourt  and  Hacker, 
1976;  Wang,  1979;  Boicourt,  1981;  Johnson  et  al.,  1983). 

Short  term  disruptions  to  routine  flow  patterns  around  the  Bay  mouth 
are  common.  Wind  forcing  of  the  shelf  waters  affects  the  nontidal  flow 
through  the  Bay  mouth.  Bottom  waters  further  than  20  km  east  of  the  Bay 
mouth  were  reported  to  flow  towards  the  Bay  with  prevailing  southerly  winds 
during  the  summer,  whereas  the  inner  shelf  surface  waters  showed  a  northward 
drift  (Boicourt  1981).  However,  strong  winds  can  produce  outflow  or  inflow 
surges.  Boicourt  (1973)  reported  a  10%  volume  reduction  in  Chesapeake  Bay 
over  48  hours  resulting  in  an  outflow  surge  extending  far  offshore. 

Nine  stations  were  monitored  for  water  quality  parameters  at  the  pro¬ 
posed  disposal  site  from  1981  through  1983  (Figure  1).  A  central  station, 

#5  (36*  59'  N  and  75*  39'  W)  and  eight  additional  stations  were  located  at 
the  cardinal  points  to  station  5:  four  stations  (16,  11,  8,  3)  were  located 
at  a  2  nm  (3.7  Km)  radius  from  the  center  and  the  remaining  four  (14,  13,  6, 
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1)  were  placed  1  nm  (1.85  Km)  beyond  the  NDS  boundary  or  5  nm  north, 
south,  east  and  west  of  center.  The  depth  of  NDS  varied  from  16  m  (52  ft) 
to  26  m  (85  ft)  . 

METHODS  AND  MATERIALS 

Water  quality  monitoring  at  NDS  stations  was  dependent  upon  ship 
availability;  however,  all  the  stations  were  sampled  for  each  season  of  the 
three  year  study  (1981  -  March,  April,  June,  August  and  October;  1982  - 
January,  April,  August,  and  October;  1983  -  February,  April,  July  and  De can- 
fa  er)  .  Duplicate  water  sanples  were  collected  at  1  m  below  the  surface  and  1 
m  above  the  bottom  in  5  or  8  1  teflon-lined  go-flo  bottles.  Aliquots  were 
withdrawn  from  the  go-flo  bottles  and  nitrite  (NO2),  nitrate  (N03),  total 
Kjeldahl  nitrogen  (TKN),  orthophosphates  (OPO^)  total  phosphorus  (TP)  chemi¬ 
cal  oxygen  demand  (COD),  pH,  turbidity,  suspended  solids,  volatile  nonfilt- 
erable  residue  (VNR),  and  chlorophylls  (_a,  j),  £,  phaeophytin  a.)  were  deter¬ 
mined. 

Field  measurements  of  temperature  (0^),  salinity  (°/oo) ,  and  dissolv¬ 
ed  oxygen  (DO)  were  monitored  at  a  surface  and  near  bottom  depth  at  each 
station  by  portable  field  meters.  The  DO  meter  was  calibrated  against  a 
Wi nkl er  titration. 

Chlorophyll  J),  £  and  phaeophytin  _a  were  measured  and  calculated  by 
the  UNESCO  method  (Strickland  and  Parsons,  1974).  Suspended  solids  and 
volatile  nonfilterable  residues  were  determined  by  drying  the  sanples  to  a 
constant  dry  weight  and  subsequently  ashing  the  samples  (APHA,  1975). 

Nitrate  levels  were  determined  by  the  Brucine  metnod  (APHA,  1975)  in 


1981  and  1982,  and  by  cadmium  reduction  (APHA,  1975)  in  1983.  Nitrites  were 
analyzed  by  sulfanilic  acid  and  (APHA,  1975).  Total  and  orthophosphate 
concentrations  were  determined  by  ammonium  molybdate  and  potassium  antimonyl 
tartrate  reactions  with  orthophosphates.  Samples  analyzed  for  total  phos¬ 
phates  were  first  digested  by  the  persulfate  method  to  oxidize  all  forms  of 
phosphorous  to  the  ortho  phosphate  form  (APHA,  1975).  Total  kjeldahl  nitro¬ 
gen  and  ammonia  samples  were  steam  distilled  and  processed  by  nessl erization 
(APHA,  1975).  Chemical  oxygen  demand  (COO)  was  performed  by  a  modified  Hach 
(1981)  procedure  in  which  samples  were  diluted  to  20  °/00  salinity  and 
treated  with  HgS04  to  reduce  chloride  interference. 

STATISTICAL  ANALYSIS 

Computer  based,  multivariate  analytical  techniques  were  employed  to 
characterize  the  spati al- temporal  characteristics  of  NDS  under  baseline 
conditions.  A  series  of  complementary  mul tivariate  statistical  approaches 
were  taken.  Discriminant  analysis  (Klecka,  1975)  was  employed  to  determine 
whether  any  significant  patterns  could  be  detected  in  the  water  quality  data 
set  with  respect  to  groups  selected  a  priori  to  the  collections.  The  groups 
were  selected  to  examine  patterns  with  depth  (e.g.,  surface  versus  bottom 
values)  and  with  season  (e.g.,  month  to  month  changes). 

In  order  to  examine  significant  patterns  detected  by  the  discriminant 
analysis  in  more  detail,  multivariate  analysis  of  variance  (MANOVA)  models 
(Hull  and  Nie,  1981)  were  employed  to  test  the  same  spatial- temporal  ef¬ 
fects.  Since  the  discriminant  models  are  often  overly  sensitive  in  terms  of 
indicating  "significant"  differences  even  when  they  do  not  occur,  the 
MANOVAs  were  considered  to  be  the  definitive  tests  of  significant  patterns 


(see  accompanying  report  entitled,  "Statistical  Significance  in  Baseline 
Monitoring"  for  a  more  detailed  discussion).  The  MANOVAs  were  also  employed 
to  test  for  any  unusual  patterns  at  NOS  under  baseline  conditions  (see  dis¬ 
cussion  of  the  season-area  interaction  model  in  accompanying  report,  ibid). 
For  each  effect,  the  MANOVAs  tested  for  significant  patterns  in  all  vari¬ 
ables  in  Doth  the  multivariate  models  and  in  univariate  analogs  to  ANOVA. 

The  third  step  in  the  analytical  procedure  was  to  employ  principal 
components  analysis  (PCA)  (Kim,  1975)  to  reduce  the  data  set  into  a  rela¬ 
tively  few  factors  which  retain  most  of  the  information  concerning  signifi¬ 
cant  patterns.  These  factors  can  then  be  plotted  for  data  presentation 
(e.g. ,  ordination) . 

Thus,  the  analytical  regime  included  a  preliminary  screening  step  (dis¬ 
criminant  analysis),  a  statistical  testing  step  (MANOVA),  and  a  data  reduc¬ 
tion/presentation  step  (PCA  ordination).  The  multiple  model  approach  to 
data  reduction,  analysis,  and  interpretation  has  been  selected  because  it  is 
recognized  that  environmental  data  sets  are  unlikely  to  conform  exactly  to 
all  of  the  assumptions  of  multivariate  statistics,  even  when  transformations 
(e.g.,  logarithmic  transformations  were  applied  to  all  data  except  pH)  have 
been  applied  in  an  attempt  to  "normalize"  the  data.  Therefore,  it  is  felt 
that  if  a  number  of  statistical  models  with  different  mathematical  approa¬ 
ches  and  assumptions  all  indicate  similar  patterns  to  be  significant,  the 
trend  can  be  considered  to  be  important.  A  second  conservative  convention 
was  adopted  in  the  interpretation  of  "significant"  trends;  only  those  rela¬ 
tionships  that  exhibited  a  high  degree  of  significance  (i.e.,  low  o  level) 
were  emphasized  so  that  the  robustness  of  the  tests  would  be  expected  to 
overcome  deviations  from  multivariate  normality.  This  convention  is  espe- 
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cially  appropriate  for  monitoring  studies,  where  numerous  statistical  tests 
may  be  run  over  time  and  Type  I  errors  (i.e.,  false  alarms)  would  be  likely 
to  occur  by  chance  alone  if  a  lower  level  of  significance  was  accepted. 
Therefore,  in  the  interpretation  of  the  findings,  only  those  trends  that  are 
highly  significant  (p  <  0.001)  or  recurrent  are  emphasized. 

In  addition  to  the  characterization  of  spatial-temporal  baseline  pat¬ 
terns,  the  water  quality  program  at  NDS  involved  the  development  of  statis¬ 
tical  techniques  for  the  estimation  of  levels  of  "minimum  detectable  im¬ 
pacts"  (MDIs):  those  levels  of  change  in  any  given  variable  which  would  be 
required  to  define  a  statistical  difference  during  trend  assessnent  studies. 
The  philosophical  concept  of  the  MDIs  and  the  methods  for  the  calculation  of 
various  MDI  models  were  detailed  in  the  accompanying  report  (ibid).  The 
MDIs  were  calculated  for  each  variable  for  single  samples,  and  for  data  sets 
(representing  data  from  a  single  post-impact  cruise)  employing  seasonal, 
baseline,  and  interaction  models. 

RESULTS  AND  DISCUSSION 

Spatio-Temporal  Patterns 

There  were  significant  changes  in  water  quality  between  each  consecu¬ 
tive  seasonal  cruise  at  the  NOS  study  area  (Table  I).  Except  for  the  Fall- 
Winter  period  of  1982-83  (October  and  February),  there  were  significant 
differences  between  the  surface  and  bottom  water  samples.  Thermal  strati¬ 
fication  was  maximum  during  the  Summer  (i.e.,  an  average  of  approximately 
10*C  difference  between  surface  and  bottom  temperatures),  and  virtually 
disappeared  during  the  colder  months  (Figure  2a).  Salinities  were  generally 
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lowest  during  the  spring  and  summer  months,  corresponding  to  periods  of 
greatest  stratification  (Figure  2b,  Table  I).  Surface  salinities  were  de¬ 
pressed  (as  low  as  25  ppt)  due  to  an  extention  of  the  Chesapeake  Bay  influ¬ 
ence  offshore.  The  effects  of  the  1980-81  drought  were  evidenced  at  the 
<N0S,  with  salinities  higher  during  the  first  year  in  comparison  to  the  re¬ 
mainder  of  the  study.  The  basic  buffering  capacity  of  the  marine  system 
held  pH  readings  near  8,  except  for  during  the  Winter  of  1982  when  values 
were  somewhat  higher  (Fig.  2c,  Table  I).  Dissolved  oxygen  levels  were 
generally  inversely  related  to  water  temperature  (Fig.  2d,  Table  I).  The 
vertical  patterns  of  oxygen  content  indicated  that  oxygen  depletion  of  bot¬ 
tom  waters  was  not  a  problem.  In  fact,  oxygen  readings  were  often  near 
saturation  levels  and  bottom  oxygen  measurements  were  only  slightly  below 
surface  values  during  two  of  the  summer  cruises. 

The  indices  of  materials  suspended  in  the  water  column  generally  exhib¬ 
ited  low  values,  but  showed  seasonal  trends  (Figure  3,  Table  I).  Lowest 
values  of  suspended  solids,  VNR  and  C.O.D.  were  generally  found  during  the 
spring  and  summer  with  highest  values  occurring  during  the  Fall.  The  major 
exception  to  this  pattern  occured  in  July  of  1983  when  higher  VNR  values 
were  observed  at  NDS,  possibly  associated  with  levels  of  dead  plant  material 
in  the  water  (as  evidenced  by  the  phaeophytin;  Figure  5d) .  The  turbidity  of 
waters  was  uniformly  low,  although  some  of  the  same  seasonal  patterns  were 
apparent  (Figure  3d) .  Highest  turbidity  readings  were  found  for  the  bottom 
waters  during  periods  of  high  chlorophyll  _a  content  (Figure  5a). 

The  nutrients  in  the  waters  of  the  study  area  were  quite  low,  often  at 
levels  below  detection  limits  (Figure  4).  Seldom  were  there  Indications  of 
significant  vertical  stratification  (Table  I).  Orthorphosphates  were  below 


Physical  parameters  monitored  at  the  proposed  Norfolk 
Disposal  Site  (NDS)  for  1981-1983.  Solid  lines  and 
asterisks  represent  surface  values,  and  broken  lines 
and  squares  are  bottom  measurements.  The  symbols 
represent  mean  values  (n*18),  with  ±  standard  error 
bars  indicated.  When  standard  error  bars  were  smaller 
than  the  symbols,  they  were0omitted  for  presentation 
purposes:  a)  temperature  (  C),  b)  salinity  (°/oo), 
c)  pH,  d)  dissolved  oxygen  (DO)  (mg/1). 
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Figure  3.  Physical  and  chemical  parameters  monitored  at  the  proposed 
Norfolk  Disposal  Site  (NDS)  for  1981-1983:  a)  suspended 
solids  (mg/1),  b)  volatile  residue  (mg/1),  c)  chemical  oxygen 
demand  (COD)  (mg/1),  d)  turbidity  of  sea  waters  (NTU). 
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detection  limits  during  1981  and  1982  and  present  in  concentrations 
approaching  the  limits  in  1983  (Figure  4a).  Total  phosphorous  was  higher 
during  the  spring  and  fall,  and  low  during  tne  summer,  although  it  remained 
below  detection  limits  during  all  of  the  1982  cruises  (Figure  4b,  Table  I). 
Likewise,  nitrites  and  nitrates  were  below  detection  limits  for  most  of  the 
study,  with  detectable  but  low  levels  reported  during  the  fall  of  1982  and 
the  spring  and  fall  of  1983  (Figure  4c,  d) .  Anmoni  a  levels  were  quite  low 
throughout  the  study,  but  a  distinct  seasonality  was  apparent:  highest 
levels  were  found  during  the  Winter  and  lowest  levels  measured  during  the 
simmer  (Figure  4c,  Table  I).  The  TKN  values  followed  a  similar,  but  consid¬ 
erably  more  depressed  pattern  with  highest  values  observed  during  Decanber 
of  1983  (Figure  4f) .  The  higher  levels  of  TKN  observed  were  believed  to  be 
associated  with  a  phytoplankton  bloom  at  the  NDS  (as  evidenced  by  the 
highest  readings  recorded  for  chlorophylls  _a  and  c_,  as  well  as,  total 
phosphorus  and  C.O.D. ;  Figures  5a,  c,  3c  and  4b,  respectively). 

Chlorophyll  content  of  the  waters  of  NDS  were  moderately  low,  but  fol¬ 
lowed  an  expected  seasonal  pattern  (Figure  5).  Chlorophylls  &  and  c  were 
highest  in  the  Fall  and  lowest  in  the  late  Spring  -  Summer  months  (Figure 
5a,  c) .  The  times  of  lowest  chlorophyll  values  corresponded  to  the  periods 
of  greatest  stratification  and  lowest  nutrient  levels.  Chlorophyll  Jo  was 
detected  in  higher  levels  in  1981  than  in  1982  or  1983,  but  no  obvious  sea¬ 
sonal  patterns  were  observed  and  the  levels  never  approached  those  observed 
for  chlorophylls  a  or  c  (Figure  5b).  Bottom  water  samples  tended  to  have 
greater  chlorophyll  levels  than  surface  samples,  except  during  the  period  of 
Spring  stratification  (Figures  5a-c) .  However,  during  the  maximum 
stratification  bottom  values  were  higher.  Phaeophytin  values  were  always 
higher  in  the  bottom  waters  with  peak  periods  of  phaeophytin  concentration 
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.  Nutrient  parameters  monitored  at  the  proposed  Norfolk  Disposal 
Site  ( MDS )  for  1981-1983:  a)  Orthophosphates  (mg/1  x  10~3), 
b)  phosphorus  (mg/1  x  10-2),  c)  Nitrite  (ug/1).  <0  Nitrate 
(mg/1  x  10~2). 
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Figure  4  (continued),  e)  Anmonia  (mg/1  x  1CT2},  f)  Kjeldhal  Nitrogen 

(mg/1  x  10"1). 
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Chlorophyll  and  phaeophytin  monitored  at  the  proposed  Norfolk 
Disposal  Site  (NDS)  for  1981-1983:  a)  chlorophyll  a^  (pg/1). 

b)  chlorophyll  (ug/1  x  10-1)»  c)  chlorophyll  c  (ug/1),  d)  Phaeo 
phytin  (ug/1 ) . 
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occurred  in  early  1982  and  mid-1983.  (Figure  5d) 

In  order  to  summarize  the  water  quality  patterns  observed  at  NDS,  a  PCA 
was  run  on  the  entire  data  set.  The  first  four  PCA  factors  accounted  for 
over  50%  of  the  variance  in  the  data  (20%,  13%,  11%,  and  9%,  respectively, 
for  factors  1-4).  The  95%  confidence  ellipses  of  the  combined  seasonal  data 
sets  were  plotted  on  graphs  of  the  principal  axes:  PCA  1  vs.  PCA  2;  PCA  1 
vs.  PCA  2;  PCA  1  vs  PCS  3;  and  PCA  1  vs.  PCA  4  (Figure  6).  The  PCA  1  repre¬ 
sented  the  changes  in  chlorophylls,  organics,  and  nutrients  that  generally 
occur  in  the  water  quality  cycle:  from  a  low  in  the  Spring  to  a  high  in  the 
Fall  and  back  again.  The  seasonal  factor  (PCA  2)  represented  the  increase 
in  temperature  and  suspended  solids,  and  the  associated  decrease  in  dis¬ 
solved  oxygen  and  certain  nutrients  which  occur  between  the  cold  and  warm 
months  (Figure  6a).  The  PCA  3  was  positively  related  to  salinity,  ammonia, 
chlorophyll  _b,  and  0.0.  which  increased  during  the  Winter;  and  negatively 
related  to  volatile  residue  which  increased  in  the  Summer  (Figure  6b). 
Primarily,  the  PCA  4  was  inversely  related  to  C.O.D.  which  was  elevated  in 
the  Fall  and  Winter;  and  to  a  lesser  extent,  PCA  4  was  directly  related  to 
chlorophyll  _b  and  to  nitrites  for  the  few  periods  when  they  were  detectable 
(Figure  6c).  Thus,  the  patterns  presented  on  Figure  6  summarize  most  of  the 
major  annual  cycles  in  water  quality  at  the  NOS,  which  have  been  detailed  in 
Figures  2-5. 

In  order  to  more  thoroughly  explore  the  spatial  patterns  of  water  qual¬ 
ity  and  to  test  trend  assessment  statistical  models  under  baseline  condi¬ 
tions,  a  series  of  MANOVAs  were  conducted  on  data  from  each  pair  of  consecu¬ 
tive  cruises  to  test  season-area  interaction  effects  (i.e.,  water  quality 
changes  occurring  within  the  NDS  not  explained  by  seasonal  patterns  and  not 
observed  elsewhere;  see  Alden,  1984  for  a  more  detailed  discussion  of  trend 


Scores  for  the  first  four  factors  of  Principal  Component 
Analysis  summarizing  the  water  quality  patterns  observed 
at  the  proposed  Norfolk  Disposal  Site  (NDS)  for  1981-1983 


assessment  statistics).  As  might  be  expected  for  baseline  conditions,  most 
of  the  tests  indicated  no  significant  differences  (a  =  0.001)  between  the 
water  quality  of  stations  within  NOS  in  comparison  with  that  observed  for 
stations  outside  the  Site.  However,  three  of  the  12  tests  (June  -  August 
1981,  August  -  October  1982,  and  October  -  February  1982-83)  indicated  that 
there  were  significant  season-area  effects. 

A  series  of  univariate  analyses  were  performed  on  the  data  from  the 
cruises  in  order  to  determine  the  parameters  responsible  for  the  signifi¬ 
cance  in  the  multivariate  models,  and  to  define  the  nature  of  the  spatial 
patterns.  The  a  priori  hypothesis  was  that  the  plume  from  the  Chesapeake 
Bay  crosses  NDS  from  time  to  time,  and  creates  spatial  differences  in  water 
quality  between  the  stations.  Normally,  the  Bay  plane  moves  in  a  pattern 
which  is  southwest  of  the  site,  but  it  has  been  observed  occasionally  to  cut 
across  this  region  (Campbell  and  Thomas,  1981).  Therefore,  a  series  of 
ANOVAs  were  conducted  to  test  whether  the  various  water  quality  parameters 
for  stations  located  in  the  southwest  portion  of  the  study  area  (1,  3,  5,  6, 
8)  were  significantly  different  from  those  observed  for  stations  to  the 
north  and  east  (11,  13,  14,  16)  (Figure  1).  A  more  definitive  test  of  spa¬ 
tial  patterns  involved  the  use  of  the  Duncan's  Multiple  Range  test  as  a 
method  to  classify  groups  of  stations  which  were  statistically  similar  with 
respect  to  each  of  the  various  water  quality  parameters.  The  univariate 
tests  were  run  on  logarithmically  transformed  data  and  Duncan's  test  was  not 
performed  unless  an  ANOVA  indicated  significant  differences  (a  *  0.001) 
between  stations  for  the  paraneter. 

In  June  of  1981,  the  stations  from  the  southwestern  portion  of  the 
study  area  had  significantly  higher  levels  of  ammonia  and  TKN  from  stations 
to  the  north  and  east  (Table  II).  The  analysis  of  August  1981  data  yielded 
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Crui se  Parameter 


S.H.  N.E. 


October  1982  Temperature 


18.37  18.64 

(0.002)  (0.003) 
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similar  results,  except  the  parameters  exhibiting  significant  elevations  at 
the  southwestern  stations  were  the  chlorophylls  and  total  phosphorous  (Table 
II).  During  October  of  1982,  the  colder  more  turbid  and  slightly  ammonia- 
enriched  waters  of  the  Bay  appeared  to  affect  the  stations  in  the  southwest¬ 
ern  portion  of  the  study  area,  particularly  in  contrast  to  the  eastern  sta¬ 
tions  which  were  influenced  by  more  oceanic  conditions.  The  February  1982 
cruise  data  provided  the  most  extensive  evidence  of  the  influence  of  Chesa¬ 
peake  Bay  water  on  the  spatial  patterns  of  the  water  quality  of  the  study 
area.  Temperature,  salinity,  dissolved  oxygen,  orthophosphates,  turbidity 
and  chlorophyll  ji  all  exhibited  patterns  which  might  be  expected  for  a  peri¬ 
od  of  plume  effects:  the  waters  of  the  southwestern  stations  had  lower  tem¬ 
peratures,  lower  salinities  and  dissolved  oxygen  content,  and  higher  turbid¬ 
ities,  phosphates  and  chlorophyll  _a,  while  the  stations  in  the  northeastern 
portion  of  the  study  area  were  more  oceanic  in  nature  (i.e.,  more  moderate 
tenperatures,  higher  salinities,  low  turbidities,  low  nutrients,  and  low 
chlorophyl  1) . 

Comparisons  to  Previous  Studies 

The  water  quality  of  the  study  area  appears  to  be  quite  good  and  gener¬ 
ally  follows  patterns  observed  by  previous  investigators.  Kester  and  Cour- 
ant  (1973)  and  Kuo  et  al.  (1975)  review  and  summarize  most  of  the  water 
quality  studies  conducted  in  the  Chesapeake  Bay  and  nearshore  coastal  wat¬ 
ers.  It  is  clear  that  the  water  quality  data  from  the  present  study  falls 
well  within  the  range  of  values  presented  in  these  reviews.  In  fact,  the 
data  reflect  the  more  neritic  location  of  the  study  area  when  compared  to 
those  reported  for  investigations  in  estuarine  regions  of  the  Chesapeake 
Bay:  nutrients,  chlorophyll  levels,  organic  content  (VNR),  suspended  solids 


and  turbidity  readings  were  lower;  salinities  and  oxygen  levels  were  higher; 
and  temperature  and  pH  were  more  stable.  As  might  be  expected,  water  qual¬ 
ity  of  the  coastal  waters  of  the  Norfolk  Disposal  Site  was  much  better  than 
that  found  for  local  estuaries  more  directly  impacted  by  man  (VSWC3,  1976). 
Likewise,  the  water  quality  conditions  in  the  vicinity  of  NDS  appeared  to  be 
much  better  than  the  coastal  regions  of  the  New  York  Bight,  with  no  indica¬ 
tions  of  the  hypoxia  or  intense  eutrophication  observed  in  those  waters 
(Pearce,  1981). 

Several  trends  were  apparent  in  comparing  the  data  from  the  present 
study  to  that  reported  by  the  Superflux  investigators  during  1980  (Campbell 
and  Thomas,  1981).  The  Superflux  program  was  conducted  during  a  record 
drought,  while  the  present  study  covered  a  period  of  much  more  normal  rain¬ 
fall  for  the  area.  It  was  not  too  surprising,  therefore,  that  watei  siality 
changes  associated  with  increased  runoff  from  the  Chesapeake  Say  were  ob¬ 
served. 

Salinities  were  generally  lower  than  those  reported  for  the  area  by 
Wong  and  Todd  (1981).  Nutrients,  particularly  ammonia,  were  observed  in 
higher  concentrations  than  those  reported  by  Wong  and  Todd  (1981).  These 
investiqators  reported  that  N/P  ratios  were  often  less  than  15,  indicating 
nitrogen  limitation  in  waters  near  the  mouth  of  the  Chesapeake  Bay.  During 
the  present  study,  a  higher  level  of  runoff  from  the  Bay  apparently  supplied 
higher  levels  of  ammonia  to  the  study  area,  so  that  phosphorous  may  be  con¬ 
sidered  to  be  more  limiting  (i.e.,  N/P  ratios  greater  than  15).  Of  course, 

the  concept  of  single  limiting  factors,  particularly  in  dynamic  coastal 
waters,  is  of  questionable  ecological  utility.  Primary  production  appears 

to  be  limited  by  a  general  paucity  of  nutrients  when  the  study  area  is  com¬ 

pared  to  the  waters  of  the  Chesapeake  Bay. 
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Suspended  solids,  volatile  residue  and  chlorophyll  concentrations  were 
generally  higher  during  the  present  study  than  during  Superflux  investiga¬ 
tions  (Gingerich  and  Oertel,  1981;  Robertson  and  Thomas,  1981;  and  Kator  and 
Zulkoff,  1981).  These  trends  also  reflect  the  higher  level  of  Bay  runoff 
during  the  period  of  the  study.  Although  still  relatively  low,  primary 
production  (as  indicated  by  chlorophyll  content)  was  presumably  higher  due 
to  the  increase  in  nutrients  from  runoff.  The  higher  levels  of  production 
may,  in  turn,  account  for  the  higher  concentrations  of  particulate  organic 
matter. 

The  horizontal  water  quality  patterns  observed  during  the  three  periods 
of  spatial  heterogeneity  (June  -  August  1981;  August  -  October  1982;  October 
-  February  1982-83)  appeared  to  echo  patterns  observed  during  the  Superflux 
study  during  June  1980  (Kendall,  1981).  The  salinity  distribution  deter¬ 
mined  by  remote  sensing  techniques  employed  during  Superflux  indicated  that 
the  Chesapeake  Bay  plume  occasionally  flowed  out  onto  the  shelf  and  cut 
across  the  area  of  the  NDS.  The  occasional  periods  of  spatial  heterogeneity 
in  water  quality  between  the  southwestern  stations  and  the  northeastern 
stations  under  more  oceanic  influences  complicate  statistical  models  to  a 
degree.  However,  the  fact  that  these  events  have  been  documented  and  de¬ 
fined  under  baseline  conditions  allows  them  to  be  taken  into  account  in 
future  trend  assessment  studies.  The  Bay  plume  effects  can  be  evaluated 
statistically  as  an  alternate  hypothesis  to  dredged  material  disposal  ef¬ 
fects,  so  that  this  potentially  confounding  natural  pattern  can  be  elimin¬ 
ated  before  the  man-made  impacts  are  evaluated. 

Minimum  Detectable  Impacts  in  Water  Quality 


In  order  to  define  the  levels  of  change  which  might  represent  "statis¬ 
tically  significant"  effects  under  the  context  of  natural  spatio-temporal 


variability,  minimum  detectable  impacts  (MDIs)  were  calculated  for  each 
variable  in  the  seasonal  data  sets  (Table  III).  The  MDIs  ranged  from  less 
than  a  5%  change  for  parameters  exhibiting  low  variability  to  over  500%  for 
the  most  variable  parameters.  However,  the  parameters  with  the  largest  MDIs 
were  generally  those  with  concentrations  very  near  to  the  detection  levels. 
Therefore,  the  absolute  values  of  the  parameters  resulting  from  the  hypo¬ 
thetical  "impacts"  were  never  extreme  in  an  ecological  sense.  Even  the 
"impacted"  values  calculated  for  single  samples  were  quite  noderate,  despite 
the  fact  that  the  level  of  change  must  be  greater  for  statistical  detection 
in  only  one  post-impact  sample. 

It  might  be  noted  that  the  direction  of  change  for  the  calculation  of 
MDIs  for  the  chlorophylls  was  chosen  to  be  negative  to  represent  inhibition 
or  toxicity  to  the  phytoplankton  population.  The  rationale  behind  this 
choice  was  that  the  chlorophyll  levels  at  NOS  were  quite  low  throughout  the 
year,  so  that  the  greatest  impact  which  could  be  reasonably  expected  would 
be  decreased  activity  by  the  primary  producers.  Unlike  the  more  estuarine 
waters  of  the  Chesapeake  Bay,  the  chlorophyll  content  (i.e.,  plant  biomass) 
of  waters  at  the  site  would  have  to  rise  orders  of  magnitude  to  even  ap¬ 
proach  levels  which  might  be  considered  eutrophic.  The  MDIs  calculated  for 
chlorophyll  increase  represented,  at  most,  a  change  of  2-3x  (200-300%). 
Thus,  although  the  chlorophyl  enhancement  MDIs  are  not  shown  in  Table  III, 
they  do  represent  a  reasonable  level  for  statistical  detection. 

The  MDIs  for  the  NDS  water  quality  indicated  that  most  changes  would  be 
statistically  detected  before  the  absolute  levels  of  the  various  parameters 
would  become  extreme  in  an  ecological  sense.  Few  of  the  projected  values 
for  "impacted"  parameters  fell  outside  of  the  natural  range  reported  by 
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Kester  and  Courant  (1973)  for  none  of  the  estuarine  of  the  lower  Chesapeake 
Bay  and  none  even  approached  the  water  quality  criteria  or  reference  levels 
recommended  by  State  and  Federal  agencies  for  the  protection  of  marine  life 
and  the  prevention  of  eutrophication  (VSWCB,  1976;  EPA,  1973).  Of  course, 
even  the  MDIs  were  somewhat  conservative  since  empirical  tests  with  data 
sets  with  computer  simulated  impacts  have  indicated  that  statistical 
detection  with  multivariate  models  is  often  possible  for  changes  of  only  30% 
of  the  estimated  MO  I  level  whenever  more  than  one  parameter  is  being 
impacted  at  the  same  time  (Alden,  1984). 

Thus,  it  is  anticipated  that  any  water  quality  impacts  associated  with 
disposal  activities  will  be  statistically  detected  during  trend  assessment 
studies  at  a  level  below  that  which  may  be  of  acute  ecological  significance. 
This  is  the  desired  situation,  if  the  monitoring  regime  and  the  associated 
statistical  models  are  to  act  as  an  "early  warning  system"  for  the  detection 
of  an  impact  before  the  environment  deteriorates  excessively  (Alden,  1984). 


Annual  Variability  and  Statistical  Detection 


As  has  been  discussed  in  a  previous  report  (Alden,  1984),  year  to  year 
variation  in  environmental  data  tends  to  confound  the  statistical  detection 
of  impacts  during  trend  assessment  studies.  The  MDIs  estimated  for  baseline 
data,  and  the  associated  variability,  are  accunul ated.  However,  as  the  data 
base  grows,  the  MDIs  would  be  expected  to  increase  asymptotically,  since  all 
but  the  most  aperiodic  baseline  "noise"  has  been  taken  into  account. 

In  order  to  determine  whether  this  theoretical  situation  occurs  with 


the  NDS  water  quality  statistical  models,  a  series  of  MDIs  were  estimated 
for  a  data  base  which  expanded  from  one  to  four  cruises  vorth  of  seasonal 
data.  The  MDIs  were  first  calculated  for  the  August  1981  data,  then  calcu- 
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lated  for  June  and  August  1981  data  sets  combined.  The  data  from  the  Summer 
cruise  of  1982  was  included  and  the  MDIs  recalculated,  and  finally,  the  1983 
Summer  data  was  incorporated  in  the  data  base  for  the  calculations. 

The  results  of  the  MDI  calculations  for  the  Summer  cruise  series  are 
presented  in  Table  IV.  Although  all  of  the  MDIs  increase  as  the  data  set 
became  more  inclusive,  the  MDIs  for  some  of  the  parameters  appear  to  have 
approached  an  asymptote:  primarily  D.O. ,  pH,  C.O.D.,  turbidity.  On  the 
other  hand,  the  MDIs  for  the  nutrients  appeared  to  still  be  increasing  with 
the  most  inclusive  data  set,  and  the  value  for  chlorophyll  a  changed  unex¬ 
pectedly  (especially  for  the  single  sample  model)  when  the  1983  data  were 
included.  When  the  series  is  expanded  to  include  the  data  now  being  col¬ 
lected  for  the  Summer  of  1984,  perhaps  more  of  the  parameters  will  have 
reached  an  asymptotic  MDI  level. 

The  trend  for  the  absolute  values  of  MDIs  to  increase  as  the  baseline 
data  set  is  expanded  may  be  disconcerting  at  first  consideration.  frie  may 
be  led  to  believe  that  it  is  better  to  monitor  baseline  conditions  for  a 
relatively  short  period,  so  that  the  MDI  levels  are  small  and  detection 
potential  is  high.  However,  this  scenario  would  lead  to  numerous  "false 
alarms"  where  natural,  albeit  somewhat  aperiodic,  changes  are  interpreted  as 
impacts.  Since  these  "false  alarms"  may  prove  to  be  costly  in  terms  of 
additional  studies  being  generated  and/or  a  moratorium  on  the  suspected 
activity,  it  would  be  more  reasonable  to  assemble  a  more  representative 
long  term  baseline  data  set.  Such  a  data  set  should  theoretically  take  into 
account  most  of  the  natural  variability  of  conditions  at  the  disposal  site, 
so  that  only  "true"  impacts  are  detected.  The  level  of  monitoring  effort 
which  produces  asymptotic  MDIs  would  be  expected  to  account  for  nearly  all 
natural  patterns,  except  those  caused  by  extreme  aperiodic  occurrences 
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(i.e.,  events  such  as  storms,  which  are  readily  documented),  or  by  long  term 
trends  (i.e.  overall  degradation  in  water  quality  in  the  region).  There¬ 
fore,  the  goal  of  baseline  monitoring  programs  should  be  to  obtain  a  suffic¬ 
ient  data  base  to  produce  asynptotic  MDIs  and  to  allow  tne  detection  of  only 

"true"  impacts  without  an  excessive  number  of  "false  alarms." 

In  order  to  demonstrate  the  effect  of  data  base  size  and  extent  on 

statistical  sensitivity,  an  example  can  be  made  of  tne  data  set  from  the 
December  1983  cruise.  Due  to  ship  time  limitations,  the  Fall  cruise  in  1983 
was  not  made  until  early  December,  a  little  over  a  month  later  than  the  Fall 
cruises  made  during  1981  and  1982.  Possibly  due  to  the  later  date  or  other 
natural  causes,  a  number  of  the  water  quality  parameters  (C.O.D. ,  TKN,  phos¬ 
phorous,  and  chlorophyll  ^)  were  at  higher  levels  than  previously  observed 
for  Fall.  When  the  October  1981  data  are  used  as  the  seasonal  baseline  in  a 
MANOVA  comparison  with  the  December  data  set,  an  extremely  high  degree  of 
statistical  difference  is  noted  (P<0.0001;  multil variate  F  =  617,  d.f  =  19, 
52).  A  comparison  of  the  canonical  variate  (c.v.)  scores  of  the  two  data 
sets  (Figure  7a)  indicated  that  all  of  the  December  values  were  outside  of 
the  99.9%  probability  limits  of  the  October  1981  set  (mean  x2  is  greater 
than  2000)  (See  Alden,  1984  for  details  of  the  probability  limits  evaluation 
procedure) . 

The  1981,  1982  and  1983  data  sets  were  analyzed  in  a  similar  manner  and 
plotted  on  the  same  c.v.  score  axis  for  comparison  purposes  (Figures  7b,  c, 
d) .  The  x2  values  for  the  99.9%  baseline  inclusion  criteria  for  the 
December  1983  data  set  dropped  dramatically  as  the  data  base  is  expanded: 
mean  x2  =  7.  52  ±  0.04  for  the  1981  baseline;  x2  ■  0.  75  ±  0.02  for  1981 
and  1982;  and  x2  *  0.043  ±  0.01  for  all  three  years  combined.  Although 
the  MANOVA  model  still  indicates  differences  (F  ■  10.21;  d.f.  »  19,472),  the 


Canonical  variate  scores  for  the  December  1983  winter  quality 
data  and  various  baseline  data  sets:  a)  scores  for  December 
1983  and  October  1981  baseline  data;  b)  scores  for  December 
1983  and  1981  annual  baseline  data;  c)  scores  for  December  1983 
and  1981-1982  composite  baseline  data;  and  d)  scores  for  December 
1983  and  1981-1983  composite  basel ine  data  (excluding  December). 


CANONICAL  VARIATE  SCORES 
4—  COO,  TKN,  TP.  Chlorophyll  fl 


probability  limits  evaluation  indicate  that  the  December  data  set  is  not 
unreasonable  for  baseline  conditions.  Since  the  use  of  the  probability 
limits  chi-square  test  is  probably  the  most  appropriate  evaluation  tool  for 
seasonal  or  baseline  model  comparisons  (Alden,  1984)  the  December  data  set 
would  not  be  considered  to  be  "impacted,"  unless  the  trend  continued  on 
subsequent  cruises. 

Therefore,  a  data  set  exhibiting  "natural"  patterns  of  water  quality 
would  have  been  considered  "impacted"  if  it  had  been  encountered  during 
trend  assessment  studies,  and  if  the  baseline  had  consisted  of  only  one 
season  or  even  one  year  of  data.  However,  when  the  data  base  is  expanded  to 
include  more  of  the  "natural"  variability  encountered  during  a  more  exten¬ 
sive  baseline  monitoring  program,  the  statistical  "false  alarm"  is  not 
sounded.  Thus,  the  importance  of  a  reasonably  extensive  data  base  is  evi¬ 
dent.  The  asymptotic  MDIs  may  provide  a  set  of  criteria  for  determining 
when  the  data  base  is  "extensive"  enough.  Obviously,  the  NDS  baseline  data 
base  is  approaching  this  level  for  certain  of  the  less  variable  parameters. 

SUMMARY  AND  CONCLUSIONS 

A  three  year  baseline  water  quality  program  was  conducted  at  the  pro¬ 
posed  Norfolk  Disposal  Site  beginning  in  1981.  The  monitoring  proqram  was 
designed  to  develop  a  data  base  against  which  future  water  quality  data  can 
be  evaluated  once  the  site  becomes  active.  The  major  goals  of  the  program 
were  threefold:  1)  to  provide  a  characterization  of  the  natural  water  qual¬ 
ity  patterns  at  NDS;  2)  to  develop  and  continually  update  multivariate  sta¬ 
tistical  models  allowing  the  detection  of  future  environmental  trends  which 
are  significantly  different  from  natural  patterns;  and  3)  to  estimate  the 
minimum  levels  of  statistically  detectable  change  in  each  parameter  (i.e. 


MDI's)  as  a  mechanism  for  the  evaluation  of  the  effectiveness  of  the  pro- 
gran. 

The  water  quality  at  the  NDS  appeared  to  be  quite  good  and  generally 
conformed  to  patterns  previously  observed.  Most  of  the  nutrients  were  quite 
low  and  many  were  below  detection  limits  throughout  most  of  the  study. 
Seasonally,  chlorophylls,  organic  load  and  most  nutrients  were  low  in  the 
spring  and  high  in  the  fall.  Suspended  solids  and  volatile  residues  were 
low  in  the  winter  and  high  in  the  sunmer,  while  dissolved  oxygen,  salinity 
and  nutrients  (particul arly  ammonia)  followed  the  opposite  pattern.  A  prin¬ 
ciple  components  analysis  was  employed  to  visually  simmarize  these  recurrent 
seasonal  cycles. 

The  overall  water  quality  was  better  than  that  reported  for  the  more 
estuarine  waters  of  Chesapeake  Bay  or  the  coastal  waters  of  the  New  York 
Bight.  Therefore,  the  region  appears  to  be  beyond  the  direct  influence  of 
many  of  man's  activities. 

The  effects  of  depth,  season  and  geographic  location  were  statistically 
defined  for  each  of  the  water  quality  parameters.  There  were  often  signifi¬ 
cant  differences  between  surface  and  bottom  measurements,  although  physical 
stratification  due  to  temperature  or  salinity  did  not  appear  to  restrict 
mixing  processes  to  the  point  of  ecological  significance.  For  example,  the 
oxygen  content  of  the  hypoliminion  did  not  drop  to  levels  that  would  prove 
detrimental  ecologically;  and  the  nutrients  seldom  showed  strong  patterns  of 
stratification.  Highly  significant  seasonal  patterns  were  defined  so  that 
they  can  be  taken  into  account  in  future  trend  assessment  statistical  mod¬ 
els.  On  the  other  hand,  differences  between  the  water  quality  at  NDS  and 
that  found  in  surrounding  waters  were  seldom  observed,  khen  significant 
patterns  were  detected,  closer  investigation  indicated  that  they  were 
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probably  caused  by  the  influence  of  the  Chesapeake  Bay  on  tne  water  quality 
of  certain  stations  (particularly  in  the  southwestern  portion  of  the  study 
area).  Like  the  seasonal  effects,  predictable  geographic  patterns  such  as 
the  pi ime-oceanic  contrasts  can  be  taken  into  account  in  the  trend  assess¬ 
ment  models. 

Minimum  detectable  impacts  (MDI's)  were  calculated  for  eacn  of  the 
water  quality  parameters.  All  of  the  MDI's  for  the  NDS  water  quality  were 
of  a  moderate  magnitude.  As  such,  statistically  significant  changes  would 
be  expected  to  be  detected  before  the  absolute  levels  would  become  environ¬ 
mentally  detrimental.  Therefore,  the  monitoring  program  in  combination  with 
the  statistical  models  would  appear  to  provide  an  effective  "early  warning 
system"  for  the  detection  of  any  impacts  which  may  be  associated  with  dis¬ 
posal  operations. 

The  effect  of  the  size  of  a  baseline  data  set  on  impact  detection  was 
explored.  As  a  baseline  is  expanded  from  one  season  to  one  year,  to  multi¬ 
ple  years,  the  MDI's  increase  rapidly  at  first  then  level  off  asynptotical- 
ly.  This  pattern  reflects  the  fact  that  more  and  more  of  the  natural  vari¬ 
ability  observed  over  time  is  being  taken  into  account  in  the  statistical 
models.  The  ultimate  implication  of  this  trend  is  that  there  should  be  a 


threshold  baseline  data  set  size  so  that  irregular  natural  patterns  are  not 
considered  to  be  "impacts."  This  effect  was  demonstrated  using  data  from 
the  NDS  monitoring  program.  The  NDS  data  base  has  revealed  the  tnreshold 
level  for  certain  of  the  parameters,  while  it  continues  to  approach  it  for 
others.  Therefore,  "false  alarms"  during  future  trend  assignment  studies 
would  be  expected  to  become  even  less  likely  as  the  data  base  is  expanded  by 
ongoing  monitoring  effects. 
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